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We derive the ground-state energy of N composite bosons made of fermion pairs using the recently 
developed composite boson many-body formalism. We concentrate on the IV-pair energy linear in 
density. We show that the scattering relevant for scattering length contains not only direct and ex¬ 
change interaction scatterings but also the dimensionless “Pauli scattering” for fermion exchange in 
the absence of fermion-fermion interaction. Numerical resolution of the resulting “ladder” integral 
equation for fermions interacting through long-range Coulomb forces — which act as effective repul¬ 
sion between excitons made of same-spin electrons and same-spin holes — shows that the prefactor 
of the IV-exciton energy linear in density is substantially decreased from its Born approximation 
value, 13-7 t/ 3, by a factor ~ 0.38 for equal electron and hole effective masses. Interestingly, this 
factor goes to zero when the hole mass goes to infinity, making the triplet-exciton gas unstable in 
this limit. 

PACS numbers: 03.75.Hh 


Although paired fermions can be approached through 
elementary fermion many-body formalisms developed in 
the 50’s, these formalisms are not suitable to handle 
fermion exchange between composite bosons made of 
fermion pairs. 

Various “bosonization” procedures have been proposed 
to transform composite bosons into elementary bosons 
interacting through effective scatterings in which a cer¬ 
tain amount of fermion exchange is included [l|. This 
can produce correct result to some problems but fails to 
address all of them, even in the dilute limit, for a very 
simple dimensional argument: while scatterings appear¬ 
ing in effective Hamiltonians are energy-like quantity by 
construction, the many-body physics of paired fermions 
is driven by dimensionless “Pauli scatterings” which orig¬ 
inate from the Pauli exclusion principle in the absence of 
fermion-fermion interaction. 

A decade and a half ago, we have proposed a con¬ 
ceptually new many-body formalism which treats com¬ 
posite bosons (“cobosons”) as entity while keeping ex¬ 
change between their fermionic components^]. This co¬ 
boson formalism relies on an operator algebra instead 
of Green function scalars. It has been mainly used to 
understand and better predict non-linear optical effects 
in semiconductors Q, these effects being controlled by 
fermion exchange in the absence of Coulomb process. 

In this Letter, we use this coboson many-body formal¬ 
ism to address the ground-state energy of N composite 
bosons made of fermion pairs, more precisely excitons 
made of same-spin electrons and same-spin holes. 

It is known that the ground-state energy of N in¬ 
teracting elementary bosons has a linear term in den¬ 
sity which comes from the repeated interaction of two 
out of N bosons. The next term in density comes 
from singular correlationjifjfl between three and more 
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FIG. 1: (a) Pauli scattering A(^P). (b) Direct Coulomb scat¬ 
tering £(™i). (c) “In” exchange Coulomb scattering 
(d) “Out” exchange Coulomb scattering £° ut (,Y). 


bosons through “bubble” processes associated with same- 
momentum transfer excitations from the condensate^]. 
In the case of composite bosons, we expect the energy 
linear in density to also come from repeated interaction 
between two out of N cobosons. The challenge is to de¬ 
termine how fermion exchange enters the associated scat¬ 
tering. 

The coboson many-body formalism generates two con¬ 
ceptually different scatterings between (i,j) and (m, n) 
states: (i) the dimensionless “Pauli scattering” A(^), 
shown in Fig. Oja), when the state change results from 
fermion exchange in the absence of fermion-fermion in¬ 
teraction; (ii) the energy-like direct interaction scatter¬ 
ing £(^), shown in Fig. [2(b), when the state change re¬ 
sults from fermion-fermion interaction in the absence of 
fermion exchange — fermions being paired in the same 
way in (i,j) and ( m,n ) states. All many-body effects 
involving cobosons can be derived in terms of these two 
elementary scatterings. 
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We here show that the scattering between (i,j) and 
(m, n) coboson pairs appearing in the “ladder” processes 
leading to the iV-coboson energy linear in density reads 
as 


— C( n i\ — cexch(nj\ _ Emn + Ejj 2-E’qq , 

Mmi/ > \mi) S \mi) ^ ’ 

(i) 

Ej,j — E 00 = E, + Ej — 2E 0 is the excitation energy of the 
coboson pair (i,j) outside the condensate, * = 0 denoting 
the single coboson ground state. £ exc/l (^) is the average 
exchange interaction scattering 


The part of 7V-coboson energy resulting from ladder 
processes between two cobosons among N can be written 
in terms of the £ scattering as 
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with £ solution of the “ladder” integral equation 


c(”° 0 ) = cO+ E 


ij^OO 


COCO 

Eqo — Etj 


(7) 


pin ( n 0\ I pout / n 3 \ 

gexch ( nj ^ S \rnij_ ' S \mi) ^ ^ j 

In interaction takes place between the “in” states 

(id) 

rd = E» ( 3 ) 

pq 

(see Fig. [2(c)) while in £°“*(^), interaction takes place 
between the “out” states (m, n); so, £°“ 4 reads as £ m with 
£ and A interchanged (see Fig. [2(d)). 

For readers knowledgeable with coboson many-body 
formalism, the scattering between (i,j) and (m, n) states 
relevant in the 7V-coboson energy cannot have a form 
different from Eq. ©■ Indeed, this scattering appears 
in matrix elements involving the fermionic Hamiltonian 
E[ between coboson states. For H acting on the right, 
the direct scattering £ leads to — £ m due to fermion ex¬ 
change when scalar products of coboson states are taken. 
By contrast, for H acting on the left, £ leads to — £ out . 
As these two exchange interaction scatterings are related 
through Pauli scattering asQ 

rc) - r\z) = (Emn Ey) Ad). (4) 

E A terms must appear along with £ m or £°“ 4 in order to 
get the same result for H acting on the right or on the 
left. Precise calculations show that they do. Then, the 
scattering given in Eq. © is the one having the required 
time reversal symmetry 

cd) = (co )*. (5) 

In the case of excitons, the effective scattering[|| pro¬ 
duced by bosonization, namely £(”■(■) — C n (mi) our 
notations, should have been rejected long ago because it 
induces a (missed) non-hermit icityQ in the Hamiltonian 
as (rO)* = £ out (/m)- Hermiticity would be easy to 
restore with £ m replaced by £ excft- ; but the EX terms of 
£ are less obvious. This once more shows the weakness 
of bosonization in problems dealing with coboson many- 
body effects. 


These two equations, along with Eq. ©. constitute the 
key results of the Letter. 

Just as the N(N — 1) prefactor in Eq. © comes from 
the number of ways to choose the two cobosons which in¬ 
teract among N , processes involving three cobosons will 
appear with a N(N — 1)(7V — 2) prefactor; and so on 
... Summing up these higher-order terms produces the 
./V-coboson correlation energy, which is expected to be 
singular as the one for elementary bosons. Its study is 
beyond the scope of the present work. 

We wish to note that, as N(N — l)/2 = 1 for N = 2 
while Eq. © is valid for arbitrary N, the part of iV-pair 
energy linear in density is related to the 2-pair energy 
through 

Sn - NE 0 ~ N{N ~ X) (S 2 - 2E 0 ) . (8) 

So, to get £n at first order in density, or equivalently the 
coboson scattering length, we can either solve the inte¬ 
gral equation © or compute the 2-pair ground state en¬ 
ergy by solving the corresponding 4-fermion Schrodinger 
equation — as previously done for short-range potential 
in the cold quantum gas context fl(il - [l2j . 

Let us now outline how the above results follow from 
the coboson many-body formalism. We look for the 
ground state of N composite bosons 

(H -£ n )\V n ) = 0. (9) 

In the dilute limit, we expect |Tv) — \&n) = Bl N \v) 
within states having some cobosons outside the conden¬ 
sate (|n) denotes the vacuum state). This leads us to look 
for | Tv) as 

I * N ) = |$JV> + E C i3 B \ B )\*N-2) + ■■■ (10) 

ij^OO 

where E>\ creates a single coboson, {H—Ei)B\ |n) = 0. To 
go further, we use two key commutators of the coboson 
many-body formalism, namely 0 

[H,B}]_ = EiBt + V?, (11) 

[V*,B]]_ = (12) 


ran 
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the direct interaction scattering £(^) being possibly re- The projection of Eq. © over |$jv), for |\I>jv) written 
placed by -£ in (”i) or even by as in Eq. (0, then reads 


= 

\mi) 


<d)-&e n d) 


(13) 


The (a, b ) indetermination at this level of calculation 
is profound. It comes from the fact that, due to fermion 
indistinguishability, there is no way to know of which 
fermion pair a composite boson is made. This unpleasant 
feature is omnipresent in the coboson many-body formal¬ 
ism because fermion exchange allows transforming a pair 
of coboson creation operators B^B^ into a sum of other 
coboson operators according to [2| 


0 = (* n \H - E n \* n ) + E (2!) 

ij^OO 


with T(° 0 i) given by 
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bIb} = -^2bIbI\(2). (i4) 

mn 

The above equation, which manifests two possible ways 
to form two cobosons out of two fermion pairs, readily 
changes £( r ”() into —£ m d) in Eq. (fT21) . Of course, phys¬ 
ical results do not depend on (a, b). 

To get H\¥ N ) in Eq. ©, we iterate Eqs. ( 11 1I12[> . This 
yields 

(H - NE 0 )\* N ) = N ^ N -V BlBl\<f> N _ 2 ) ~C° 0 ) ■ 

mn 

(15) 

and, for (i,j) ± (0,0) 

(H - E^ - (TV - 2)E^B\B]\<S> N - 2 ) 

( 16 ) 

pq 

within terms in (N — 2) and higher which correspond 
to processes involving more than two cobosons excited 
from the condensate. These higher-order processes lead 
to correlation energy, as previously explained. So, we 
shall not consider them. 

Next, we project the Schrodinger equation © onto 
|$jv) and B^BW^n- 2 )■ These projections make ap¬ 
pear scalar products of coboson states. Starting from 
= N\Fn [see Ref. 0], we get 

(®N- 2 \B n B m \$ N ) — N\F N _ 2 [s m 0 S n0 — A(£°)) , (17) 


for A n =£ n —NE 0 . To go further, we note that 

E (<w<w - a(^))h(«)= az) - rd) (23) 

pq 


whatever (a, b) taken for H('()), due to the fact that two 
successive fermion exchanges between two cobosons re¬ 
duce to an identity. 

Equations (1151171) give the Hamiltonian mean value 

(■ H) n as 


($ n \H\$ n ) 

{®n\®n) 


~ NEn + 


N(N - 1) F n _ 2 
2 F n 


(cO-r©) 


(24) 

in agreement with previous work[13|. We also find from 

Eq. ©ED 


£ n ~ (H) n + ^ E (C(oD + AjvAC'))^,- (25) 

N 00 

where, with the help of Eq. ©, the C(mi) scattering is 
just the one defined in Eq. ©. 

To obtain Cjj, we project Eq. © on B i B i\$ N - 2 ). 
Using Eqs. (fl6l) and (fT51) . we find, after some algebra, 

(E 00 - E mn + A N )c mn ~ N{N ~ 1} (c^^+AivAdo)) 

+ E (Cd) + AjvAd)) Cij (26) 

b# oo 


and 

(^jv-2 \B n B m B\ | <l>7v—2) 

— (A — 2 )\F n . 2 ( 5 mi 5 nj - A(^) + (* <—► jj) , (18) 

within terms in (N — 2). These scalar products follow 
from the iteration of the other two key commutators of 
the coboson many-body formalism, namely 

[B m ,B^]_ = 5mi - Dmi , (19) 

[-Dmi , B]\ _ i = E B l (Ad) + (* <-► 3)) • (20) 


Next, we drop the Ajv terms: indeed, the change A n in 
the energy of N cobosons induced by interaction depends 
linearly on coboson number or equivalently on sample 
volume. So, the Ajv terms provide contributions to £n 
with improper volume dependence. They come from dis¬ 
connected processes which cancel out, as standard in in¬ 
teraction expansion. 

By setting c^- = dijN(N — l)/2, the above equation 
reduces to 

Cdo) = (Eoo-E mn )d mn ~ Cdo)+ E Cd)dii > (27) 

bV oo 












As F n _ 2 /F n ~ 1 at first order in density 0, this gives 
the 7V-coboson energy as 
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£n ^ NE 0 + N{N 2 1 } C@ (l + -A) (28) 

where T denotes the correction factor to the Born ap¬ 
proximation 


T = 



Y C(o! iKj 
oo 


1 \ ' C(oi)C(io) 

C(oo) E 00 - E ij 


(29) 


Equation (JHD readily follows from the above equations. 

In the following, we concentrate on Coulomb potential 
to address the ground-state energy of N excitons made of 
same-spin electrons and same-spin holes. As previously 
shown[2|, the direct Coulomb scattering £(°o) is equal to 
zero because repulsion between two electrons or between 
two holes is as strong as attraction between one electron 
and one hole. By contrast, the exchange Coulomb scat¬ 
tering £ m (oo)i equal to £ 0 “*(oo) due to Eq. (g}, is given 
by (ob/L) D ^_d in 3 D Rydberg units, with £3 = — 267 r /3 
and £2 = —( 87 t — 3157r'V512) for electrons and holes in 
the same quantum well|14|. C(oo) thus reduces in 3 D to 
-f n ® = ( 267 T /3 )(a B /L) 3 . 

This analytical result has already been obtained by 
Keldysh and KozlovEEJ in the late 60’s. One great 
advantage of the coboson many-body formalism over 
Bogoliubov-like procedure used up to now[l 6 j, is to pro¬ 
vide a physical understanding to this term: it comes from 
one exchange interaction inside the condensate with di¬ 
rect scattering reducing to zero in the case of Coulomb 
potential. To the best of our knowledge, explicit result for 
higher-order Coulomb processes has not been derived yet. 
According to Eq. CO. these processes contain direct and 
exchange Coulomb scatterings as well as Pauli scatter¬ 
ing multiplied by the excitation energy of the pair states 
at hand. Here also, the coboson many-body formalism 
allows catching the physics of these Coulomb processes 
readily. 

To calculate contribution coming from more than one 
Coulomb interaction, we first note that the exciton pair 
(i,j) excited from the condensate, which enters the re¬ 
peated interaction of Eq. ED, has a total center-of- 
mass momentum equal to zero, Qj + Q_, = 0 , because 
Coulomb and exchange processes conserve momentum. 
These pairs can have arbitrary relative motion indices 
(yi^Vj). However, due to energy denominators, domi¬ 
nant terms come from excitons staying in the relative 
motion ground state, vq. So, the scatterings CiZii) f° r 
the dominant processes in the ladder equation ED reduce 

to cr§:;^"§*^) = c(Q m ,Qi). 

We have numerically solved the integral equation ED 
with C(m») replaced by C(Qm, Qi)- This amounts to only 
keeping the Is ground state vq for the exciton relative mo¬ 
tion while the center-of-mass momentum is treated using 



FIG. 2: Correction factor T to the Born approximation as a 
function of electron-to-hole effective-mass ratio, m e /vrih- 


a dense mesh in Q. Details of the numerical implementa¬ 
tion is similar to what is described in Ref. El Numerical 
results for the correction factor T to the Born approxi¬ 
mation, defined in Eq. (l29l) . are presented in Fig. [2] as a 
function of electron-to-hole effective-mass ratio ( m e /mh )• 
It is seen that T monotonically decreases from T ~ —0.62 
for m e /rrih = 1 down to F = — 1 when m e /rrih — > 0, indi¬ 
cating a complete cancellation of the density-linear term 
in this limit. 

This cancellation remains true when all relative motion 
levels v are kept in the (i,j) sum of Eq. ED- This can 
be shown analytically by considering Eq. E3 for m = 
(Q, ^ 0 ) and n = (—Q, isq). The excited pair energy E mn — 
Eqo goes to zero when the hole mass, ie., the center-of- 
mass mass, goes to infinity. So, in this limit, the (i,j) 
sum in Eq. (l27l) must go to — £( q'0“ 0 ) which, for Q->0, 
reduces to — C(8o)- According to Eq. (l29l) . F thus goes to 
— 1 in the large hole mass limit. This physically means 
that, when the hole mass approaches infinity, the exciton 
gas does not resemble the |\H) state given in Eq. (TTOl) but 
should be better described as a set of localized excitons. 
This limit will be studied elsewhere. 

To conclude, we have derived the energy of N com¬ 
posite bosons linear in density at all orders in fermion- 
fermion interaction, using the coboson many-body for¬ 
malism which treats the Pauli exclusion principle in an 
exact way. This linear term follows from a “ladder” in¬ 
tegral equation that is numerically solved in the case 
of Coulomb potential between same-spin electrons and 
same-spin holes making semiconductor (triplet) excitons. 
We show that higher-order interactions drastically re¬ 
duces the value of the density-linear term obtained within 
the Born approximation, down to zero in the infinite hole 
mass limit. This coboson many-body formalism opens 
the route to securely address the far more complex next- 
order term in density expected to behave as n 3//2 instead 
of n 2 [see Ref. EH- This formalism can also be used for 
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other than Coulomb potential and for other fields than 
Semiconductor Physics, such as Cold Quantum Gases. 
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